function clean_EEG = f_artifact(raw_EEG,m,fs)
% first, set the initial 5 sec to be 0. (opening artifact)
raw_EEG(:,1:round(5*fs)) = 0;

B =  jadeR(raw_EEG,m);
clean_EEG.data = B*raw_EEG;
clean_EEG.A = raw_EEG/clean_EEG.data;
